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Collective mass tensor derived from the cranking approximation to the adiabatic time-dependent 
Hartree-Fock-Bogoliubov (ATDHFB) approach is compared with that obtained in the Gaussian 
Overlap Approximation (GOA) to the generator coordinate method. Illustrative calculations are 
carried out for one-dimensional quadrupole fission pathways in 5 Fm. It is shown that the collective 
mass exhibits strong variations with the quadrupole collective coordinate. These variations are 
related to the changes in the intrinsic shell structure. The differences between collective inertia 
obtained in cranking and perturbative cranking approximations to ATDHFB, and within GOA, are 
discussed. 
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I. INTRODUCTION 

Microscopic understanding of nuclear collective dy- 
namics is a long-term goal of low-energy nuclear theory. 
Large amplitude collective motion (LACM), as seen in 
fission and fusion, provides a particularly important chal- 
lenge. Those phenomena can be understood in terms of 
many-body tunneling involving the mixing of mean fields 
with different symmetries. We have yet to obtain a mi- 
croscopic understanding of LACM that is comparable to 
what we have for ground states, excited states, and re- 
sponse functions. 

For heavy, complex nuclei, the theoretical tool of 
choice is the self-consistent nuclear density functional 
theory (DFT) 0, Q. The advantage of DFT is that, 
while treating the nucleus as a many-body system of 
fermions, it provides an avenue for identifying the es- 
sential collective degrees of freedom and provides an ex- 
cellent starting point for time-dependent extensions. The 
time-dependent Hartree-Fock-Bogoliubov (TDHFB) the- 
ory appears, in principle, to provide a proper theoretical 
framework to describe the LACM. However, the main 
drawback of TDHFB, when applied to fission, is its in- 
ability to describe the quantum-mechanical motion under 
the collective barrier. 

On the other hand, the adiabatic approximation to 
TDHFB (ATDHFB) has been successfully applied to the 
LACM J3H12J]. The main assumption behind ATDHFB, 
well fulfilled in the context of spontaneous fission, is that 
the collective motion of the system is slow compared to 
the single-particle motion of individual nucleons |l|, 1 131 ) . 
According to the path formulation of the fission prob- 
lem 14j , ATDHFB provides the best framework to tackle 
the problem of nuclear dynamics under the barrier. An- 
other advantage of ATDHFB is that it provides a con- 
nection between the microscopic many-body theory and 
phenomenological models based on collective shape vari- 



ables. 

The main theoretical input for an estimate of fission 
half-lives is collective inertia (mass tensor) and collective 
potential. ATDHFB provides the best framework to cal- 
culate mass tensor [14]]. However, in most applications, 
various approximations are adopted. In the commonly 
used cranking expression, for instance, the derivatives 
with respect to collective coordinates (i.e., collective mo- 
menta) are evaluated using the perturbation theory, and 
the Thouless-Valatin self-consistent terms yielding time- 
odd fields are neglected. The resulting collective masses 
are known to be too small [6, 9]; hence it is imperative 
to go beyond the perturbative cranking treatment. 

In the self-consistent investigations of Ref. Q , based on 
the Gogny energy density functional, collective masses 
were calculated by explicitly evaluating the collective- 
coordinate derivatives appearing in the ATDHFB mass 
expression. The resulting collective mass obtained in 
such an approach turned out to exhibit appreciable vari- 
ations along the collective path, suppressed in the per- 
turbative cranking treatment. Furthermore, they noted 
that ATDHFB cranking mass could be an order of mag- 
nitude greater than the perturbative cranking mass. As 
noted in Ref. [9j, the enhanced masses obtained in the 
improved analysis can significantly impact the calculated 
fission lifetimes. 

The main goal of this work is to investigate the AT- 
DHFB cranking mass using the nuclear DFT approach 
with Skyrme energy functionals. The paper is organized 
as follows. Section [H] summarizes the basic ATDHFB 
expressions for collective inertia obtained in Ref. [6j. 
The approximate cranking, perturbative cranking, and 
Gaussian Overlap Approximation (GOA) formulations 
are given in Sec. IIIII The illustrative examples of cal- 
culations are contained in Sec. IIV1 where the results are 
presented for 256 Fm. Finally, the main results are given 
in Sec. IVl 



II. ATDHFB THEORY 

This section contains a brief derivation of the collec- 
tive mass tensor in the ATDHFB framework. Although 
some of the expressions are well documented in the lit- 
erature [l|, Q , we repeat them here for the sake of com- 
pleteness with particular attention paid to various ap- 
proximations involved. 



In Eq. (j6|) and in the following, the quantities in the 
canonical basis are denoted by symbols with breve ac- 
cents [ia |. 

The HFB energy matrix E in the canonical basis is 
non-diagonal and is given by 



Env — Cj, (h - A) M „ - rj A M p s* 



(9) 



where 



A. Summary of HFB 

We begin with the HFB approach. In what follows, we 
use the same notation as in Ref. [15[ . The HFB formalism 
can be conveniently expressed in terms of the generalized 
density matrix, 1Z, defined as 



n 




K 
1-p* 



(1) 



where p and k are the particle and pairing densities and 
1Z 2 = 1Z. The energy variation results in the HFB equa- 
tion 

[w,n] = o , (2) 

which can be written as a non-linear eigenvalue problem: 



W 



where 



A B* 
B A* 



W = 



A B* 

B A* 



h-X A 

-A* -h* + X 



E 
-E 



(3) 



(4) 



E is a diagonal matrix of quasiparticle energies E^, X 
is the chemical potential, and matrices h and A are the 
particle-hole and pairing mean- field potentials [1] , respec- 
tively. 

For the sake of comparison with the commonly used 
BCS formalism, it is quite useful to write the HFB equa- 
tions in the canonical representation. The single-particle 
canonical wave function |/i) can be expanded in the orig- 
inal single-particle (harmonic oscillator) basis \n) as 



Im> = Y1 Dn » I' 1 )' 



(5) 



where the unitary transformation D is obtained by diag- 
onalizing the density matrix p. In the canonical basis, 
the HFB wave function is given in a BCS-like form: 



-™-lii> — ^fiVfii/ , t> ^y — S ^V^O^^, 



U,. = Un 



Vn 



(6) 

(7) 



where the phase s^, for the time-even quasiparticle 
vacuum considered here, is defined through the time- 
inversion of the single-particle states 



?» = -vIm) 



(8) 



f]% = u n v v ± u v v,j, and £±„ = u^Uv =p v^v v . (10) 

The diagonal matrix elements of the matrix E^ u can be 
written as [H, [l5| : 



E„ 



( V - A) 2 + A 2 ^ 



(11) 



Even though the above equation resembles the BCS ex- 
pression for quasiparticle energy, it involves ft M/i and A^ , 
which are respectively obtained by transforming the HFB 
particle-hole and the pairing fields to the canonical basis 
via the transformation ([5]). It is only in the BCS ap- 
proximation that these quantities can be associated with 
single-particle energies and the pairing gap. 



B. Summary of ATDHFB 

The ATDHFB approach is an approximation to the 
time-dependent HFB theory, wherein it is assumed that 
the collective velocity of the system is small compared to 
the average single-particle velocity of the nucleons. The 
generalized HFB density matrix is expanded around the 
quasi-stationary HFB solution 1Z$ up to quadratic terms 
in the collective momentum: 



K = TZ Q + K x + K 2 , 



(12) 



with IZi being time-odd and TZq and 1Z 2 time-even den- 
sities. The corresponding expansion for the HFB Hamil- 
tonian reads 



w = w + yvi + w 2 



(13) 



Employing the density expansion (fT2")l . the HFB energy 
can be separated into the collective kinetic and the po- 
tential parts. In terms of the density expansion ([12")) , the 
kinetic energy is given by 

K = ^Tr(W n 2 ) + ^TiiW 1 n 1 ) 

= Jli(*o[72o,Wi])-i([Wa,^o][Wo,Wo])<14) 

In the usual ATDHFB treatment, the second term in- 
volving 1Z 2 is neglected, and the kinetic energy can be 
written in the familiar form: 



IC=^q 2 M, 



(15) 



where the collective mass is given by 



M = ^Tr(ft [fto,fti] 



2q \ dq 



(16) 

(17) 



The trace in the above expression can easily be evaluated 
in the quasiparticle basis. To this end, one can utilize the 
ATDHFB equation " 



iTl = [Wo,1li} + [W l ,Tl } 



(18) 



In the quasiparticle basis, the matrices TZq , Wo , Wi , 1Z\ , 
and IZq are represented by the matrices Q, £q, £ i, Z, and 
IF, respectively: 



TZo = AGA\ 

Wo = A£ a A^ , 
Wi = A8iA\ 
K x = AZA^ , 
7t = AF./^ , 



(19) 

(20) 
(21) 
(22) 
(23) 



where 



A 



A B* 

B A* 



(24) 



is the matrix of the Bogolyubov transformation, and 



Q = 




1 



£o = 



E 
-E 



(25) 



ATDHFB equation (TT8")) can now be written in the quasi- 
particle basis as 



iF=[£ ,Z} + [£ 1 ,g] 



(26) 



This 2x2 matrix equation is, in fact, equivalent [6J to 
the following equation, 



iF = EZ + ZE + Ei 



(27) 



where the antisymmetric matrices F , Z , and E\ are re- 
lated to J 1 ", Z, and £\: 



T 



F 

-F* 



z 
-Z* 



[£i,S\ 



Si 
-£? 



(28) 



(29) 



In the case of several collective coordinates {ft}, the 
ATDHFB equation ([18]) must be solved for each coordi- 
nate. 






and the collective mass tensor becomes: 



i ( dlZn , 4. 



(30) 



(31) 



Then, in terms of the corresponding matrices F % and Z^ , 
the collective mass tensor is given by 



M 



-Tr| F"Z' J -F l Z r 



2ft ij V 



(32) 



The expression (|3"2")) for the mass tensor contains the ma- 
trix Z % , which is associated with time-odd density matrix 
1Z\ and can, in principle, be obtained by solving the HFB 
equations with time-odd fields. The time-odd fields have 
been incorporated in mass-tensor calculations only in a 
limited number of cases. For instance, in Ref. [ljj, time- 
odd fields have been included in the HF study with a 
constraint of cylindrical symmetry. The time-odd fields 
have also been incorporated in the HFB study in an ap- 
proximate iterative scheme with the collective path based 
on the Woods-Saxon potential [g|. 



III. APPROXIMATIONS TO ATDHFB 

This section contains the summary of various com- 
monly used approximations to the exact ATDHFB ex- 
pression ([52"]) . 

A. Cranking approximation 

In most of the studies, the time-odd interaction matrix 
Ei appearing in Eq. (|27| is neglected. In the following, 
this approximation will be referred to as the cranking 
approximation (ATDHFB-C). In the absence of the term 
involving E\ , the Z-matrix can be easily obtained in the 
quasiparticle basis from the equation: 

-iFJ u , = (E lt + E v )Z i l „, (33) 

and the collective cranking mass tensor is given by: 

i (pi* pj i pi pj*\ 



M) 



2qiQj 



En + E V 



(34) 



It should be noted that Eq. (|33|) is diagonal in the quasi- 
particle basis and not in the canonical basis. The es- 
sential input to the ATDHFB-C mass tensor (13~4"|) is the 
matrix F. In the following, F is evaluated in both canon- 
ical and quasiparticle basis. 



1. Canonical basis 



To begin with, Eq. (|23p can be written explicitly in 
terms of the HFB eigenvectors: 



fto 



d 



Po 




(35) 



AFB T - B*F*A j < AFA T - B*F*B^ 
BFB T - A*F*A* BFA T - A*F*B^ 



Evaluating the matrix elements of (|35|) in the canonical 
basis, we obtain 



one arrives at an approximate "BCS-equivalent" expres- 
sion for the matrix elements of F in the canonical basis: 



«;/ 



{UpVv + V^U V ) 



dpo 



(36) 



)IU 



By differentiating the HFB equation [Wo , TS-o] = with 
respect to %, the derivative of the density matrix in 
(|55| can be expressed in terms of the derivatives of the 
particle-hole and the pairing mean-fields. The resulting 
2x2 matrix equation is 



a%^ao 



£o,J- 



= 0. 



(37) 



By employing the properties of h and A with respect to 
time reversal, we obtain 



dh* 

dqi 
/<9A* 



V dq l 



\IV 



fiv 



S ^v% 



dh 



dA 
dqi 



fiiy 



fj,y 



(38) 
(39) 



and by approximating the HFB energy matrix in the 
canonical basis by its diagonal matrix elements, 



E„ 



Ofij/J^fx i 



(40) 



F, 



\IV 



Efx + E v 



«^(/»'-aV + ^(aV » (4i) 



where x 1 = dx/dqi with x = h, A, or A. In the following, 
the results obtained by using this approximation will be 
called ATDHFB-C C . 

Using relations (|4"0)l and (J4*T|) . the collective mass ten- 
sor (|34p can now be expressed in terms of the derivatives 
of the mean-field potentials with respect to the collective 
coordinates qi and BCS-like quasiparticle energies (fTTj) . 



M 



c° 



2q i q i ^ 



En* pj | pi T?j* 
flV (IV ' fJ,V (IV 

E^ + Ey 



(42) 



In the one-dimensional case, the resulting expression 
agrees with that of Ref . [9] . 



2. Quasiparticle basis 

In order to obtain the expression for matrix F in the 
quasiparticle basis, we invert Eq. (|23p and write the ma- 
trix expression for T = A^IZqA 



T = 



A^pqA + A^k B - B^k* A - B^p* B A^p B* + A^k A* - B^k* B* - fit p *A* 
B T Po A + B T k B - A T k* A - A T p* B B T Po B* + B T k A* - A T k* B* - A T p* A* 

I 



(43) 



Elements (1,1) and (2,2) of J" vanish because TZq is pro- 
jective, IZq = TZq. Equating the above expression with 
Eq. ([25]). we obtain 



The simplex structure of various terms in Eq. (J44I) is 
given by 



- F* = B 1 p A + B 1 k B - A 1 k* Q A - A 1 p* B . (44) 

In the following, we evaluate the above expression in the 
simplex basis, as the mean-field analysis has been per- 
formed by imposing this symmetry. In this basis, the 
HFB wave function has the following structure 



B 



B+ 
B_ 



and A 



A+ 
A- 



(45) 



The density matrices acquire the following forms in the 
simplex basis 



B* + Bl 




B*B T 



P+ o 
p_ 



B\A T _ \ _ ( k+ 
B*Al J ~ { k_ 



(46) 

(47) 



B T PQ A = 
A T PoB = 
B T k B = 



This yields: 




B T po-A_ 



A T + P* 0+ B + 




A T k* A 



Blpo+A+ 


A T _p* _B_ 



B + k 0+ B- 


A T _k* _A+ 



~ F* = 



F+ 
F_ 



(48) 



(49) 



where 



F_ 



F + = Bip 0+ A+ - Aip* _B. 



BiKiQ-^-B— — A_ 



A 



+ ■ 



Bipo-A- - A\pl + B^ 
B^k -B + - Alko + A. 



(50) 
(51) 



Since F is antisymmetric, we have obviously F+ = — F-, 
which is fulfilled explicitly provided k+ = — k_. 



3. Calculation of derivatives 



neglecting the time-odd interaction terms in the ATD- 
HFB equation and off-diagonal matrix elements of the 
HFB energy matrix (j4"0)) , the derivatives are not eval- 
uated explicitly but are obtained using a perturbativc 
approach. A complete description of the perturbativc 
cranking model as applied to the nuclear fission process 
can be found in Refs. [l7l - l2T| . 

The perturbative cranking expression for the mass ten- 
sor is obtained by approximating the mean-field deriva- 
tives in Eq. (j3"6")l by canonical-basis expressions. For in- 
stance, the matrix element of h l can be approximated 
by 



The collective mass involves either derivatives of the 
density matrices or the mean-field potentials. It should 
be stressed that these derivatives must be calculated in 
the original single-particle basis \n) as the canonical ba- 
sis §5$ varies with {qi\. In the following, we show how to 
evaluate the collective derivatives in the one-dimensional 
case of single collective coordinate, the quadrupole de- 
formation q. To this end, we approximate the derivative 
of the density operator p or k at a deformation point 
q = qo by means of the Lagrange three-point formula for 
unequally spaced points qo — Sq, qo, and qo+Sq' [9l. [l3 . [l6j : 



dp 



-p( qo -6q) + -^-J-p(q ) 



Sq(Sq + Sq') 
Sq 



5qSq' 



-prp{qo + 5q'). 



(52) 



8q'{dq + Sq')' 

The reason for the use of unequally spaced grid in 
Eq. (|52[) is purely numerical: the constrained HFB equa- 
tions cannot be precisely solved at a requested deforma- 
tion point q. 

The corresponding matrix element in the canonical ba- 
sis can be expressed through the matrices D nu of the 
canonical transformation <j5j> : 

dp\ —Sq' 

Oq 



ia> 



Sq(Sq + Sq' 

Sq ~ Sq' 2 
SqSq' V » ^ 
Sq 



7T Yl D n lfi (p(<l0 - Sq)) niri2 D n2V 



(53) 



5q'(Sq + Sq 



77 Yl D n^(p(lo + Sq')) ni n 2 D n2V . 



It should be noted that the canonical matrix D nv in the 
above expression corresponds to the deformation point, 
q , at which the mass is evaluated. Furthermore, as men- 
tioned above, the density matrices at the three deforma- 
tion points in (|52")l need to be calculated using the single- 
particle basis \n) with the same basis deformation. 



B. Perturbative cranking approximation 

The perturbative cranking approximation (ATDHFB- 
C p ) has been widely used for the evaluation of the col- 
lective mass tensor. In this approximation, apart from 



(i/|h'V)«('»/»-MW/'>> 
for /i =/= v which, together with 



(54) 

(55) 
(56) 



for h,. 



/i MM and A M 



expression for the cranking mass tensor 
(mIWH^'Im) 



-A M pS^, leads to the following 



AC~£ 



\>* 



t(h u -K)*{E u + E„) 



(v) 



o i 



2£„ 



where 



(57) 



I^'i — en 

v- ~ w ~ 2E 2 



[&M - A*) - {K - A)AjJ. (58) 



Assuming a weak state dependence of A^ [22J , neglecting 
the derivatives of A and A [21j , and using the identity 



v 

Vv 



(hp - A)A„ - (h v - A)A M 



E^K 



E„A^ 



(59) 



one arrives at the following perturbative cranking mass 
tensor: 



A<«£ 



(p\W\v)(v\V\p) 
{E u + E v f 



+ ^2 



(r? + 



(60) 



where the sums run over the whole set of canonical states. 
This expression resembles the standard cranking expres- 
sion for the collective mass tensor [13, LL8|, 120, uLLf origi- 
nally derived for a phenomenological mean field h. 



C. Gaussian overlap approximation 

To compare cranking expressions with those obtained 
within the GOA, it is convenient to introduce the S ma- 
trices 231: 



s£° = 



y {p\hy){ V \hl\p) ^ + 



jl.ll 



(E„ + E„) 



K 



^iv t 



(61) 



It is immediately seen that for the mass tensor of Eq. (|60j) 
one has M c = S^\ In the case of GOA, also assum- 
ing weak state dependence of pairing and neglecting the 



derivatives of A and A M ^ , one obtains [23J, |2J 



(62) 



M GOA = S^ [s(i)l £( 2 



Evaluating the matrix elements of h l entering Eq. (|6"Tj) 
perturbatively, one can express S explicitly through the 
matrix elements of the constraining field operators Qi. 



S (K) 



Afd 



M ( 



K) 



Afd 



(63) 



where the energy-weighted moments M^ K ' are given by 



M 



(K) 






in' 



{E„ + £„ 



(64) 



D. Treatment of proton and neutron contributions 

The above expressions for the mass tensor are valid 
for one kind of fermions only. In the case of the cranking 
approximation, the total mass tensor is a sum of neutron 
and proton contributions: 



-Hotal = M 



C 



■Mi 



(65) 



In the GOA, however, the total inverse inertia 
(A-^^j ) _1 for a composite system is given as a sum of 
proton and neutron inverse covariant inertia tensors [24J : 



(^Ci)- 1 



(M 



GOA\-l 



{M; 



GOAwl 



(66) 



IV. RESULTS 



The illustrative calculations were performed for the 
nucleus 256 Fm by using the SkM* energy density func- 
tional 25] in the particle-hole (ph) channel. In the 
particle-particle (pp) channel we employed the density- 
dependent pairing interaction in the mixed variant of 
Refs. SHI: 



V T {f) = V TQ {l-p{f)/2p Q )8{i f ) 



(67) 



where r = n,p and po = 0.16 fm _1 . To test the accuracy 
of various approximations, we carried out both HF+BCS 
and HFB calculations. The pairing interaction strengths, 
which were adjusted to reproduce the neutron and proton 
ground-state pairing gaps in 252 Fm, are (in MeV fm 3 ): 



V n0 



^372.0 , V, 



p0 



-438.0. 



(68) 



The fission pathways were studied in the previous 
Ref. [H[ using the SkM*-HF-BCS approach with the se- 
niority pairing interaction. It has been found that the 
SkM* energy density functional favors the asymmetric 



fission pathway in 256 Fm, and our HFB results are con- 
sistent with this result. The one-dimensional collective 
pathway, determined by the axial quadrupole moment 
Q = Q20, was obtained by means of the HFB solver 
HFODD [21. 
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FIG. 1: (Color online) The quadrupole mass parameter (top: 
total; middle: neutron contribution; bottom: proton con- 
tribution) along the fission pathway of 256 Fm calculated in 
SkM*+HF+BCS as a function of the mass quadrupole mo- 
ment. The ATDBCS-C results (triangles) are compared with 
those obtained in the perturbative cranking approximation 
(ATDBCS-C P , circles) and Gaussian overlap approximation 
(ATDBCS GOA , diamonds). See text for details. 

Using the self-consistent solutions along the fission 
pathway, we calculate the collective quadrupole mass 
parameter using various approximations described in 
Sec. IIIII First, we discuss results obtained within the 
HF+BCS formalism (dubbed ATDBCS). Figure [T] com- 
pares the results of the non-perturbative cranking ap- 
proach (ATDBCS-C) with the perturbative cranking ap- 
proximation (ATDBCS-C P ) and Gaussian overlap ap- 
proximation (ATDBCS GOA ). In ATDBCS-C the deriva- 
tives of the density matrices and the mean- field potentials 
have been obtained using the Lagrange formula, which 
requires the knowledge of self-consistent solutions in sev- 
eral neighboring deformation points. We have evaluated 
the density matrices for quadrupole deformations ranging 
from Qiq= to 320 b in steps of 1 b. The derivatives were 
obtained by using the 3-point Lagrange formula (1531) , and 



also the 5-point Lagrange formula [l(j. The results for 
collective mass obtained with 3-point and 5-point expres- 
sions differ only in the third decimal place; hence, in the 
following, we shall stick to the 3-point Lagrange formula. 
It needs to be stressed that - in order to guarantee consis- 
tent labeling of canonical states - the underlying single- 
particle basis should be identical for all three points in 
Eq. (|53[) involved in the derivative evaluation. This has 
been achieved by performing HF+BCS calculations using 
the same basis deformation for all neighboring points. 

As seen in Fig. Q] the total ATDBCS-C mass exhibits 
a rather irregular behavior characterized by the presence 
of several sharp maxima. Some of these peak-like struc- 
tures, although considerably suppressed, also show up in 
ATDBCS-CP and ATDBCS GOA . 
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FIG. 2: (Color online) Total HF+BCS energy E tot (top), pair- 
ing energies E pp (total, proton, and neutron; middle), and the 
HF energy E p h (bottom) calculated along the fission pathway 
for 256 Fm. The borders between different self-consistent con- 
figurations are marked by vertical solid lines. The dotted lines 
mark positions of peaks in the collective ATDBCS-C mass pa- 
rameter of Fig. [T] 

To unravel the origin of the peak structures in the 
collective mass, the total energy of 256 Fm is depicted 
in Fig. [21 together with corresponding pairing-energy 
(Epp) and HF energy (E p h) contributions. The one- 
dimensional total energy curve shows several discontinu- 
ities due to intersections of close-lying energy sheets (sur- 
faces) with very different mean fields. The correspond- 
ing pathways can in fact be well separated when studied 



in more than one dimension of the collective manifold 
|28| . The diabatic jumps between various energy sheets 
have been disregarded when computing the collective in- 
ertia shown in Fig. [1] Indeed, in such cases the adiabatic 
theory is unable to provide a meaningful result for the 
collective mass. 

There are also configuration changes within each path- 
way. Because of pairing correlations, these changes are 
adiabatic in character [13,|3(|. Still, they manifest them- 
selves in the collective inertias through the appearance of 
peaks [3ll433| . Figures Q] and [2] nicely illustrate this point: 
the peaks in the collective ATDBCS-C mass parameter 
appear in the regions of large local variations in E pp and 
-Eph that are indicative of changes in the shell structure 
with elongation. (We note that the local variations in 
the total energy are much weaker than those in pairing 
and HF energies, due to the well-known anticorrelation 
between pairing and HF energies, clearly seen in Fig. [5]) 

Two general conclusions can be drawn from the results 
of Fig. [U First, the ATDBCS-C P and ATDBCS GOA in- 
ertia show fairly similar behavior, with the ATDBCS- 
CP mass being systematically larger. Second, the exact 
treatment of derivative terms in ATDBCS-C gives rise 
to less adiabatic behavior in the corresponding collective 
mass. 

The results obtained within the HFB framework are 
presented in Fig. [31 The results obtained in the 
canonical approximation ATDHFB-C C (03]), perturba- 
tive treatment of derivatives ATDHFB-C P ^7^, and 
ATDHFB GOA (|62p were obtained by using the canonical 
HFB wave functions and employing the diagonal ( "equiv- 
alent BCS") ansatz. The ATDHFB-C calculations (JM]) 
were carried out in the full quasiparticle basis. 

The most interesting finding is that the collective 
mass in ATDHFB-C is very close to that obtained 
in ATDHFB-C C . Similar to the HF-BCS case, the 
ATDHFB-C p and ATDHFB GOA results follow each other 
with the ATDHFB-C P mass being systematically larger. 
Again, the exact treatment of derivatives gives rise to less 
adiabatic behavior of collective mass that manifests itself 
through the presence of peaks. 

In Ref. 0, the quadrupole collective mass was evalu- 
ated in the canonical basis and exhibited a singular be- 
havior at certain deformation points. The primary reason 
for this singularity is due to the pairing collapse at certain 
deformations that results in unphysical phase transition 
and the presence of unavoided level crossings. In our 
work, the peak structures are present at nonzero pairing 
and are related to the shell structure changes along the 
fission pathways. 

In order to highlight the differences between HFB and 
BCS treatments, Fig.0]shows the quadrupole masses ob- 
tained in these approaches in the region of the ground- 
state minimum and the inner fission barrier of 256 Fm 
{Q20 £ [20, 100]). This region plays a crucial role in the 
evaluation of fission half-lives. It is evident from Fig. @] 
that the non-perturbative cranking masses ATDBCS-C 
and ATDHFB-C have very similar behavior. On the 



, 256 Fm 



Total 




40 80 120 160 200 240 280 

Quadrupole moment Q 20 (b) 

FIG. 3: (Color online) The quadrupole mass parameter 
(top: total; middle: neutron contribution; bottom: pro- 
ton contribution) along the fission pathway of 256 Fm calcu- 
lated in SkM*+HFB as a function of the mass quadrupole 
moment. The ATDHFB-C results (triangles) are com- 
pared with those obtained in the canonical approximation 
(ATDHFB-C C , squares), perturbative cranking approxima- 
tion (ATDHFB-C P , dots), and Gaussian overlap approxima- 
tion (ATDHFB GOA , diamonds). See text for details. 



other hand, the masses calculated in the perturbative ap- 
proximations, ATDBCS-C p and ATDHFB-C P , are quite 
different for Q20 < 40. Furthermore, in this region, the 
perturbative masses appear to be quite large as compared 
to the cranking values. 

The high-frequency fluctuations of collective mass can 
be traced back to the imperfect numerical convergence 
of HFB calculations. In the present work, we assumed 
the accuracy of 0.001 MeV for the total energy. This re- 
sults in an uncertainty of about 0.002 h 2 /(MeV b 2 ) in 
the collective inertia. If required, the precision of these 
calculations can be increased at the expense of an appre- 
ciably higher CPU time. 



V. SUMMARY 

The primary motivation of the present work has been 
to assess various approximations to the collective mass 
for fission. The collective mass plays a crucial role in 
determining the adiabatic collective motion of the nu- 
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FIG. 4: (Color online) Similar as in Fig. [3] except for 
ATDHFB-C (filled triangles), ATDBCS-C (open triangles), 
ATDHFB-C P (dots), and ATDBCS-C P (circles) in the nar- 
rower region of 20 b <Q20<100b. 



cleus and strongly impacts predicted half-lives. In the 
majority of previous studies, cranking approximation to 
collective mass has been employed, in which the time- 
odd fields are ignored and the collective momenta (i.e., 
derivatives with respect to collective coordinates) needed 
in the evaluation of the ATDHFB mass are calculated us- 
ing the perturbation theory. 

In our study, we performed the full ATDHFB cranking 
treatment of quadrupole inertia. The numerical evalua- 
tion of the derivatives appearing in ATDHFB mass ex- 
pression poses a serious computational challenge as the 
accurate self-consistent HFB solutions need to be ob- 
tained at several neighboring points around every de- 
formation along the fission pathway. By comparing 
three- and five-point approximations, we conclude that 
the three-point Lagrange formula provides a reasonable 
description of collective derivatives. 

The main conclusions of this work can be summarized 
as follows. 

• The collective masses obtained in non-perturbative 
treatment of derivatives show more variations due 
to shell structure changes along the fission path as 
compared to the perturbative approximation and 
GOA. 

• The collective mass in full ATDHFB-C is very close 



to that obtained in ATDHFB-C C and ATDBCS-C. 
This means that the diagonal approximation P0|) 
for the HFB energy matrix is a very reasonable one. 

• The ATDHFB-CP and ATDHFB GOA inertias ex- 
hibit very similar pattern, with the ATDHFB-C P 
mass being systematically larger. A similar conclu- 
sion has been reached for the HF+BCS case. 

• The main difference between HFB and HF+BCS 
calculations shows up in the perturbative treat- 
ment: the collective masses calculated in ATDBCS- 
C p and ATDHFB-C P are sometimes fairly different. 

• Considering the differences between exact cranking 
results and ATDHFB-C P and ATDHFB GOA vari- 
ants, we conclude that the perturbative treatment 
of derivatives cannot be justified. 

The present work deals with the cranking approxima- 
tion to ATDHFB in which only time-even mean fields 
have been kept when evaluating the collective inertia. 
The discussion of the full ATDHFB treatment, including 
the time-odd response that is expected to play a signif- 



icant role in the description of collective dynamics [34 
will be the subject of a forthcoming study. 



Acknowledgments 

This work was supported in part by the National Nu- 
clear Security Administration by the National Nuclear 
Security Administration under the Stewardship Science 
Academic Alliances program through DOE Grant DE- 
FG52-09NA2946f; by the U.S. Department of Energy 
under Contract Nos. DE-FG02-96ER40963 (University 
of Tennessee), and DE-FC02-09ER41583 (UNEDF Sci- 
DAC Collaboration); by the NEUP grant DE-AC07- 
05ID14517 (sub award 00091100); by the Polish Ministry 
of Science under Contracts Nos. N N202 328234 and 
N202 231137; and by the Academy of Finland and Uni- 
versity of Jyvaskyla within the FIDIPRO programme. 
Computational resources were provided by the National 
Center for Computational Sciences at Oak Ridge Na- 
tional Laboratory and the National Energy Research Sci- 
entific Computing 



[1 
[2 

[3; 
[4: 

[5 
[6. 

[r\ 

[8 
[9 

[io; 
in 

[12 

[13 

[14 

[is; 

[16. 

[17; 
[is; 

[19 



Vautherin, Lecture 



Quentin, and M. 



P. Ring and P. Schuck, The Nuclear Many-Body Problem 
( Springer- Verlag, Berlin, 1980). 

M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. 
Phys. 75, 121 (2003). 

M. Baranger and M. Veneroni, Ann. Phys. 114, 123 
(1978) . 

S.J. Krieger and K. Goeke, Nucl. Phys. A 234, 269 
(1974). 

D.M. Brink, M.J. Giannoni, and M. Veneroni, Nucl. 
Phys. A 258, 237 (1976). 

J. Dobaczewski and J. Skalski, Nucl. Phys. A 369, 123 
(1981). 

F. Griimmer, K. Goeke, and P.-G. Reinhardt, Lecture 
Notes in Physics 171, 323 (1982). 
P.-H. Heenen, H. Flocard, and D. 
Notes in Physics 171, 338 (1982). 
E.Kh. Yuldashbaeva, J. Libert, P. 
Girod, Phys. Lett. B 461, 1 (1999). 

J. Libert, M. Girod, and J.-P. Delaroche, Phys. Rev. C 
60, 054301 (1999). 

M.J. Giannoni and P. Quentin, Phys. Rev. C 21 
(1980). 

M.J. Giannoni and P. Quentin, Phys. Rev. C 21 
(1980). 

W. Nazarewicz, Nucl. Phys. A 557, 489c (1993). 
J. Skalski, Phys. Rev. C 77, 064610 (2008). 
J. Dobaczewski, W. Nazarewicz, T.R. Werner. 
Berger, C.R. Chinn, and J. Decharge, Phys. Rev. 
2809 (1996). 

M. Abramowitz and LA. Stegun, Handbook of Mathemat- 
ical Functions (Dover, New York, 1970). 
S.T. Belyaev, Mat. Fys. Medd. Dan. Vid 
11) (1959). 

D.R. Bes and Z. Szymariski, Nucl. Phys. 28, 42 (1961). 
S.G. Nilsson, C.-F. Tsang, A. Sobiczewski, Z. Szymariski, 



[20] 
[21] 

[22] 
[23] 

[24] 

[25] 

[26] 



2060 


[27 


2076 


[28] 




[29] 


J.F. 




C 53, 


[30 




[31 


,emat- 


|32 




|33 


L (No. 


[34 



S. Wycech, C. Gustafson, I.-L. Lamm, P. Moller, and B. 
Nilsson, Nucl. Phys. A 131, 1 (1969). 

A. Sobiczewski, Z. Szymahski, S. Wycech, S. G. Nilsson, 
J. R. Nix, C. F. Tsang, C. Gustafson, P. Moller, and 

B. Nilsson, Nucl. Phys. A 131, 67 (1969). 

M. Brack, J. Damgard, A.S. Jensen, H.C. Pauli, V.M. 
Strutinsky and C. Y. Wong, Rev. Mod. Phys. 44, 320 
(1972) . 

K. Sieja and A.Baran, Phys. Rev. C 68, 044308 (2003). 
A. Staszczak, S. Pilat, and K. Pomorski, Nucl. Phys. A 
504, 589 (1989). 

A. Gozdz, K. Pomorski, M. Brack, and W. Werner, Nucl. 
Phys. A 442, 26 (1985). 

J. Bartel, P. Quentin, M. Brack, C. Guet, and H.B. 
Hakansson, Nucl. Phys. A 386, 79 (1982). 
J. Dobaczewski, W. Nazarewicz, and M.V. Stoitsov, Pro- 
ceedings of the NATO Advanced Research Workshop 
The Nuclear Many-Body Problem 2001, Brijuni, Croa- 
tia, June 2-5, 2001, eds. W. Nazarewicz and D. Vretenar 
(Kluwer, Dordrecht, 2002), p. 181. 

J. Dobaczewski, W. Nazarewicz, and M.V. Stoitsov, Eur. 
Phys. J. A 15, 21 (2002). 

A. Staszczak, A. Baran, J. Dobaczewski, and W. 
Nazarewicz, Phys. Rev. C 80, 014309 (2009). 
J. Dobaczewski et at, Comput. Phys. Commun. 102, 
166 (1997); 102, 183 (1997); 131, 164 (2000); 158, 158 
(2004); 167, 214 (2005); 180, 2361 (2009). 
J.W. Negele, Nucl. Phys. A 502, 371c (1989). 
J.J. Griffin, Nucl. Phys. A 170, 395 (1971). 
G Schiitte, Z. Phys. A 296, 331 (1980). 
P. Arve and GF. Bertsch, Phys. Lett. B 215, 1 (1988). 
K. Matsuyanagi, M. Matsuo, T. Nakatsukasa, N. Hino- 
hara, and K. Sato, J. Phys. G: Nucl. Part. Phys. 37, 
064018 (2010). 



